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ABSTRACT 

Using relativistic, axisymmetric, ideal MHD, we examine the outflow from a disk around a compact object, 
taking into account the baryonic matter, the electron-positron/photon fluid, and the large-scale electromagnetic 
field. Focussing on the parameter regime appropriate to 7-ray burst outflows, we demonstrate, through exact self- 
similar solutions, that the thermal force (which dominates the initial acceleration) and the Lorentz force (which 
dominates further out and contributes most of the acceleration) can convert up to ~ 50% of the initial total energy 
into asymptotic baryon kinetic energy. We examine how baryon loading and magnetic collimation affect the 
structure of the flow, including the regime where emission due to internal shocks could take place. 

Subject headings: gamma rays: bursts — ISM: jets and outflows — MHD — methods: analytical — relativity 



1. INTRODUCTION 

Cosmological gamma-ray bursts (GRBs) evidently involve 
the emission of the isotropic equivalent of ~ 10 51 - 10 54 ergs 
on a typical time scale of a few seconds in the vicinity of a 
newly formed stellar-mass black hole or rapidly rotating neu- 
tron star (see reviews by Piran 1999 and Meszaros, Rees, & 
Wijers 1999). To avoid strong attenuation of the 7-rays by pair- 
production interactions with lower-energy photons (the "com- 
pactness" problem), the Lorentz factor 7 of the emitting gas 
must exceed ~ 10 2 - 10 3 . The most natural way of powering 
the burst seems to be the conversion of the kinetic energy of the 
relativistic outflow into nonthermal (most likely synchrotron) 
radiation. Because of the high initial energy density and char- 
acteristic photon energy, the dynamical evolution of GRBs has 
traditionally been modeled in terms of the expansion of an ini- 
tially opaque electron-positron (e ± ) fireball. Although a pure- 
ed 1 fireball could in principle attain 7 > 10 3 in the optically thin 
region where the observed emission must originate, the fraction 
of the energy carried by pairs would be only ~ 10~ 5 (Grimsrud 
& Wasserman 1998). If, however, the fireball contains enough 
baryons (of total mass Mb) that it becomes "matter dominated" 
before it gets optically thin, then most of its initial (suffix z) en- 
ergy £f is converted into bulk kinetic energy of baryons (with 
asymptotic Lorentz factors 7^ w £j/Mt,c 2 ) instead of escap- 
ing as radiation. Baryon loading of the fireball at the source 
is expected in view of the highly super-Eddington luminosity 
involved, and viable models in fact face the issue of why it is 
not so efficient as to render 700 too low (the "baryon contami- 
nation" problem). 

According to currently accepted GRB formation scenarios, 
the burst is powered by the extraction of rotational energy from 
the central black hole or neutron star, or, alternatively, from the 
debris disk left behind when the mass near the origin collapses 
into a black hole. In either case, strong (> 10 14 G) magnetic 
fields provide the most plausible means of extracting the energy 
£ i on the burst time scale. If the field were initially weaker, 
it would likely be rapidly amplified by differential rotation or 
dynamo action to the requisite strength. The magnetic energy 
might be dissipated near the origin in a series of flares, giv- 
ing rise to a "magnetic" fireball (e.g., Narayan, Paczyriski, & 
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Piran 1992; Meszaros, Laguna, & Rees 1993). Alternatively, 
the magnetic field may have a large-scale, ordered component 
that could help guide and collimate the outflow, and, if it is 
strong enough, also contribute to its acceleration (e.g., Usov 
1994; Thompson 1994; Meszaros & Rees 1997; Katz 1997). 
Even if the flow is not magnetically driven, the field may be 
strong enough to account for the observed synchrotron emis- 
sion (Spruit, Daigne, & Drenkhahn 2001). Magnetic collima- 
tion is, in fact, consistent with the observational indications for 
GRB jets (e.g., Sari, Piran, & Halpern 1999) and could be very 
helpful in reducing the source energy requirements to plausible 
values. Furthermore, if the outflow is largely Poynting flux- 
dominated, then the implied lower radiative luminosity near the 
origin could alleviate the baryon contamination problem. Mag- 
netic fields have thus come to be regarded as the favored means 
of driving GRB outflows. 

The purpose of this Letter is to provide an overview of the 
magnetohydrodynamics (MHD) of GRB outflows, clarifying 
the relationship between the thermal (fireball) and Lorentz ac- 
celerations, identifying the parameter regimes where qualita- 
tively different behaviors are expected, and demonstrating that 
the observationally inferred properties of these outflows can be 
attributed to magnetic driving. A detailed exposition of the 
semianalytic solutions on which this overview is based is given 
in N. Vlahakis & A. Konigl, in preparation (hereafter VK). 

2. THE MHD DESCRIPTION 

For the purpose of illustration, we assume that the outflow 
originates in a debris disk around a stellar-mass black hole (e.g., 
Meszaros & Rees 1997; Fryer, Woosley, & Hartmann 1999). In 
light of the considerations of § 1 , we concentrate on the case of a 
magnetically driven, axial jet that has a sufficiently high baryon 
loading to insure that it is matter dominated when it becomes 
optically thin (see §4), and in which a significant fraction of 
the Poynting flux is eventually converted into kinetic energy. 
As discussed by Spruit et al. (2001), such configurations can 
be validly described by the MHD approximation on the spatial 
scales of interest, and, furthermore, they are not expected to dis- 
sipate a substantial amount of magnetic energy along the way. 
Correspondingly, we adopt the equations of relativistic, ideal 
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MHD to describe the flow. We anticipate that, near the origin, 
the thermal energy associated with the radiation and pairs 
is nonnegligible and that the optical depth is large enough to 
ensure local thermodynamic equilibrium. We therefore assume 
that the gas (consisting of baryons with their neutralizing elec- 
trons as well as of photons and pairs) evolves adiabatically with 
a poly tropic index of 4/3. Assuming a quasi-steady poloidal 
magnetic flux function A and changing variables from (A , I , t ) 
to (A , I , s = ct — £), with t the arclength along a poloidal field- 
line, it can be shown that all terms with derivatives with respect 
to s are negligible when the flow is highly relativistic. As is 
elaborated on in VK, the equations are then effectively time in- 
dependent and the motion can be described as a frozen pulse 
whose internal profile is specified through the variable s (see 
also Piran 1999). We discuss time-dependent effects in §6. 

Assuming also axisymmetry, the full set of MHD equations 
can be partially integrated to yield several field-line constants: 
the mass-to-magnetic flux ratio ^>a(A,s), the field angular ve- 
locity il(A , s) (which equals the matter angular velocity at the 
footpoint of the field line in the disk), the total (kinetic + mag- 
netic) specific angular momentum L(A , s), the total energy- 
to-mass flux ratio p{A,s)c 2 , and the adiabat Q(A ,s) = P/p^ 
(where po and P are, respectively, the rest-mass density and to- 
tal gas pressure measured in the fluid frame). Two integrals 
remain to be performed, involving the Bernoulli and trans- 
field force-balance equations. There are correspondingly two 
unknown functions, which we choose for convenience to be 
x = 05£l/c, the radial distance of the field line [in cylindrical 
coordinates (z,G3,0)] in units of the "light cylinder" radius, 
and the "Alfvenic" Mach number M = \/4npo~£(~/V p / B p ) = 
\/£/47rpo (where £ is the enthalpy-to-rest energy ratio, and 
where V p and B p are the poloidal components of the velocity 
and magnetic field, respectively, measured in the central ob- 
ject's frame). The solutions derived in VK are obtained under 
the most general ansatz for radial self-similarity [in spherical 
coordinates (r , 9 , <j>)], in which the shape r(9,A) of the poloidal 
field lines is given as a product of a function of A times a func- 
tion of 6: r = Ti{A)T 2 {0) (or, equivalently, A = A[m/G(9)]; 
Vlahakis & Tsinganos 1998). 

Initially, the outflow evolves as a fireball: it expands isotropi- 
cally in the comoving frame, with the large-scale field only act- 
ing as a guide (e.g., Meszaros et al. 1993). If the baryon loading 
is large enough, this continues until all the thermal energy of 
the radiation field and the pairs is transformed into kinetic en- 
ergy of the baryons. Subsequently, the Lorentz force converts 
Poynting flux into baryon kinetic energy, and it also collimates 
the flow. The Lorentz force remains dominant (implying that 
the field is force free) all the way to the point where the ki- 
netic energy of the baryons becomes comparable to the electro- 
magnetic energy. We assume that the asymptotic shape of the 
flow is cylindrical and that the final stage is characterized by an 
equipartition between the Poynting and kinetic energy fluxes. 
(The latter assumption is motivated by the exact "cold" solu- 
tions of Li, Chiueh, & Begelman 1992.) The transition from a 
parabolic to a cylindrical shape occurs far downstream from the 
Alfven point (9 <C 1), where the ratio of the kinetic energy flux 
£7 2 poc 3 to the Poynting flux cEB^/An (where E is the elec- 
tric field amplitude) equals M 2 /x 2 . If $ is the (small) opening 
half-angle of the outflow, then | VA | sin(6>-i9) = OA/ dr. Sub- 
stituting this into the Bernoulli equation, we find the slope of 
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where a = A£l 2 /(c 3 ^) is the magnetization parameter. 
Equipartition between the Poynting and kinetic energy fluxes 
(M 2 = x 2 ) in the final cylindrical stage {dXS = 0) means that 
rdA{r ,9)/Adr = p/{2a). Substituting this into equation (1) 
and specializing to the force-free (M 2 <C x 2 ) parabolic regime, 
gives <f Inz/c/lnOS = 2, so the shape of the force-free field lines is 
z/G5j = (05/05,) 2 /(2tanz9,), with the initial opening half-angle 
(cf. Contopoulos 1995). 

3. SCALING LAWS 

Beyond the light surface (which is close to the Alfven sur- 
face, since the flow is nearly force free), the large-scale electro- 
magnetic field is given by 

(B Z7 B m ,B < p 7 E Z7 E m7 E ( i ) ) = - - A 2^> - l> 2z' -1 '^) ' ^ 

The comoving field amplitude satisfies B 2 Q = B 2 -E 2 (a Lorentz 
invariant) and B 2 a «Bjt^/7 2 « B 2 ^, which together imply 

Bcf, w E m (or, equivalently, LVL rs pc 2 ). 

We approximate the outflow as a pair of shells that move 
in opposite directions from the debris disk, each having an 
initial meridional cross section (Az),- x (A05)„ with (A05), = 
10 6 (Aa3) (i 6 cm comparable to 05, = 10 6 C5 (i 6 cm (the mean 
radius of the debris disk) and (Az), w cAt = 3 x 10 11 Ar : 
cm (where At = lOAfi s is the total duration of the burst). 
The total baryonic mass of outflowing matter is Mj, = 
4 7 rEi,(Aa3),(Az) i 7, / 9oi = 2 x 10- 7 7,A)/ ! 2ro, !6 (Aa3), !6 Af 1 M Q , 
where po; = 10 2 poi,2 g cm" 3 . The total energy is = pMbC 2 = 
4 x 10 51 (M/ 1 4 )7iP0(,2roi.6(AC3) !i 6Afi ergs, and initially it re- 
sides predominantly in the electromagnetic field; the initial 
thermal energy (associated with the enthalpy of the photons 
and pairs) is &M/,c 2 = (£ ( //u)£,-. 

The constancy of the mass-to-magnetic flux ratio = 
A-Kpo^Vp/Bp and equation (2) imply 

7Po03 2 = 7;/°o;G5 2 . (3) 
If is the temperature in units of the electron rest energy, 
then, by the conservation of specific entropy along the flow, 
P/(po<d) = const, or (since P oc 4 ) 

p = POi (e/e i ) 3 . (4) 

Energy flux conservation, in turn, gives 

^-(mnB^/v A c 2 ) = p. (5) 

So long as £ > 1 (the fireball phase), £ oc 0. Also, by equa- 
tion (5), £7 is constant since the specific Poynting flux does not 
vary along the flow for a force-free field. Equations (3) - (5) 
thus imply 

7 = 7,(03/03,), po = po,(03,'/03) 3 , = 0,(05, /05) . (6) 
These scalings are the same as in a spherically symmetric, hy- 
drodynamic (HD) fireball, except that the cylindrical radius 03 
replaces the spherical radius r in the magnetically guided case. 

Beyond the point where £ = 1, magnetic acceleration takes 
over and the Lorentz factor continues to increase until the cylin- 
drical regime (03 = 05oc) is reached. The magnetic acceleration 
is associated with a small (linear in 03) deviation of the Poynting 
flux from its force-free value, 
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By combining equation (7) with equations (5) and (3), it is seen 
that the scalings of 7 and po with 03 given by equation (6) are 
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maintained also in the magnetic acceleration zone. In fact, since 
the density in the central object's frame is 7/50 = 7,po;(OJ,/G3) 2 
and the 05-width of the shell (given the parabolic shape of the 
field lines) evolves as A05/(A0J), = 03/03,, mass conservation 
implies that, while the flow continues to expand isotropically 
in the comoving frame, the width of the shell remains constant 
[Az = (Az),] in the central source's frame (consistent with the 
frozen-pulse approximation). 

4. BARYON LOADING 

Two parameters determine the basic character of the flow: 
The first is the initial enthalpy of the flow in units of the rest 
energy of the baryonic component, = 400 6f (po;,2)~\ where 
we approximate the initial radiation field by a blackbody distri- 
bution and the e ± pairs by a Maxwellian distribution. 
The second is the initial Poynting flux in units of the baryon 
rest- energy flux. As we focus on Poynting flux-dominated out- 
flows, this ratio is (j,— w fi. 

Initially the flow is optically thick. The distance where it be- 
comes optically thin is determined by The pertinent regimes 
of this parameter were previously analyzed for the case of a ra- 
dially expanding HD fireball (see Piran 1999). We now adapt 
this analysis to the present (nonradial) geometry. 

The opacity has two contributions: one due to pairs (r±) and 
another due to the electrons that neutralize the baryons (rfc). As 
the flow expands, its temperature decreases. When 6 w 0.04, 
the pair number density becomes negligible and the opacity 
from there on is determined by the neutralizing electrons. For 
the radiation-dominated case, tj, < 1 at that point, correspond- 
ing to £• > &± (with e± = 2.5 x 10 7 92(AC5k 6 ). If 6 < 6±, the 
flow continues to be optically thick until Tb = 1 . If at that point 
the radiation energy density exceeds the baryon rest-energy 
density, then one is still in the radiation-dominated regime. This 
corresponds to & > with £ a = 6 x 10 3 6* /3 (AG3)^ 3 . For 

< £a,, the flow is matter dominated when it becomes opti- 
cally thin. Finally, £, w 1 correspond to the cold limit. Since 
magnetic acceleration now plays a role, the outflow is not nec- 
essarily nonrelativistic in this limit: only if fi w 1 also holds will 
the flow be Newtonian. The most plausible regime for GRB 
outflows is 1 < < £ib, and we adopt it from here on. 

5. STRUCTURE OF THE FLOW 

As the flow moves outward it passes through the following 
points (see Fig. 1): 

1. 03 1 = 03, : origin of the outflow, 7,- w 1. 

2. 032 = 03,- y/3/2: slow-magnetosonic surface, V p w c/y/3. 

3. 033 = c/fi : Alfven (w light) surface. 

4. 03 4 = (c/fl)y^/l 



classical fast-magnetosonic surface, 
•yVp s=s y / (B 2 — E 2 )/(4'KpQ^). Unlike the purely radial case 
(Michel 1969), this point is located at a finite distance from 
the origin. The bulk of the magnetic acceleration occurs down- 
stream from this point through the "magnetic nozzle" mech- 
anism (Li et al. 1992). In essence, this term is a shorthand 
for the fact that an MHD flow can continue to accelerate un- 
til it crosses the modified fast-magnetosonic singular surface. 2 
This effect is not purely relativistic: it is manifested in an ex- 
act solution of the nonrelativistic MHD equations where all the 
singular surfaces (including the modified-fast one) are crossed 
(Vlahakis et al. 2000), and there, too, most of the acceleration 

2 The modified (and not the classical) fast-magnetosonic surface is a singular 
and transfield equations (e.g., Bogovalov 1997). 



occurs downstream of the classical fast-magnetosonic surface 
(but upstream of the modified-fast surface). 

5. 03 5 =256,03, : r± = 1. 

6. 036 = £;03, : end of thermal acceleration, 7 w The entire 
initial thermal energy of the photons and pairs has by this point 
been converted into baryon kinetic energy. Magnetic accelera- 
tion effectively starts here. 

7. 03 7 = 5 x 10 3 (A03)-f pj^GS/ : r = n = 1, the flow becomes 
optically thin to Compton scattering. 

8. 038 = 03oo w (p/2)03,-: cylindrical flow regime. Near- 
equipartition between the Poynting and the baryon kinetic- 
energy fluxes is attained, with 7 = 700 w /i/2. 



10* |r 

10 3 

10 2 
10' 
10° 
10" 

10 12 

I 

N 

10 10 
10 8 
10 6 



-(atfffi/F A c) 



(a) 



A 



y 



-••••••>-" 



B 

II 

B 



(b) 



10' ra(cm) 10° 



10° 



10" 



FIG. 1. — Exact self-similar solution of the relativistic MHD equations, il- 
lustrating the GRB outflow structure discussed in this Letter, (a) The Lorentz 
factor 7, the ratio £ of the enthalpy to the rest energy, and the ratio of the 
Poynting flux to the rest-energy flux (top curve) are shown as functions of 03, 
the distance from the axis of rotation, along the innermost field line, (b) The 
meridional projections of the innermost and outermost field lines are shown on 
a logarithmic scale, along with a sketch of the black-hole/debris-disk system. 
The vertical lines mark the positions of the various transition points listed in 
§5 along the innermost field line. 

The inequalities 03i < 032 < 033 < 034 < 05g always hold, 
as do, in the regime of interest, OS5 < OS7 (valid for < 
£,±) and 036 < 037 (valid for < £,7,). The validity condi- 
tion for the inequality 037 < 03g, which expresses the require- 
ment that the flow be optically thin to Compton scattering 
in its final (cylindrical) stage, corresponds to an upper bound 

on the baryon loading: po/,2 < (m/10 4 ) (A03)^ 6 [or po/,2 < 

0.4(A03)7i03:g /3 (Af 1 )- 2 / 3 (£:,/10 51 ergs) 2 / 3 ]. We note, however, 
that the increase of 7 with 03 is slower than linear as the cylin- 
drical regime is approached (see Fig. 1), which implies that 
03 oo is larger than (p/2)03,- and the actual upper bound on the 
initial baryon mass density is larger than the analytic estimate 
(typically by a factor of a few). 

All of the above points are indicated in Figure 1, where 
an illustrative solution of the steady, relativistic MHD equa- 
tions is presented. The figure demonstrates the validity of the 

surface for the steady MHD equations when one solves simultaneously the Bernoulli 



4 



scaling 7 w 05/05, over several decades in 03 and the separa- 
tion of the thermal (7 < £,) and magnetic (7 > £,) acceleration 
regimes. It also manifests the significant collimation from the 
comparatively large initial opening half-angle (i?,- w 24°) to a 
very nearly cylindrical geometry (attained on scales > 10 14 cm 
from the origin). Equation (7) for the Poynting flux is also ver- 
ified. It is, furthermore, seen that approximate equipartition 
(7 s=s -05flB ^ / ac 2 ) holds during the final phase of the flow. 

6. APPLICATION TO GRBs 

To demonstrate how the model outflows discussed in §5 
could account for observed GRBs, we adopt as fiducial source 
parameters a 7-ray fluence of ~ 10" 5 ergs s" 1 and a distance 
<~ 3 Gpc, which imply an isotropic equivalent energy of <~ 
10 53 ergs. If the 7-ray emitting material is in fact confined 
to a pair of cones of opening half-angle ^3°, then the actual 
radiated energy is ~ 1.5 x 10 50 ergs. As we noted in §1, the 
emission is believed to be powered by the kinetic energy of the 
relativistic outflow, and the most common interpretation is in 
terms of internal shocks produced by the collision of overtaking 
shells (see Piran 1999). Assuming that the kinetic-to-radiative 
energy conversion efficiency is ~ 10% and that the magnetic 
stresses can transfer as much as <~ 1/2 of the initial energy into 
baryon kinetic energy, we infer £, « 3 x 10 51 ergs. If the en- 
ergy is deposited in the form of a Poynting flux over Af w 10 s, 
then the field at the origin must be ~ (^c/oj^AOJXAffi 2 ) 1 / 2 w 
3 x 10 14 G for 03* w (AGJ)i w 10 6 cm and Q w 10 4 s" 1 (repre- 
sentative of a debris disk that extends to the last stable orbit 
of a rotating solar-mass black hole). To satisfy the constraint 
®7 ^5 03s (see §5), the maximum baryonic rest-mass density is 
Pm,2 ~ 1. In the solution presented in Figure 1, we adopted as a 
plausible temperature 0, w 1/ \pl, corresponding to w 100. 

The variability properties of the observed GRBs have been 
interpreted in terms of N = 100 N2 distinct shells that are ejected 
with slightly different Lorentz factors and whose subsequent 
collisions give rise to the "internal" shocks responsible for the 
7-ray emission (e.g., Piran 1999). Adopting 7 oc 05 oc y/z, 
one may integrate the equation of motion for each shell and 
show that two neighboring shells starting with A7; ~ 1 will 



not collide for as long as they move in the flow acceleration 
zone. Only after the shells reach the constant-velocity, cylin- 
drical flow regime, will the two shells collide (at a height 
z f w 7^(Az),/7V = 7.5 x 10 lf V/10 4 ) 2 Afi//V 2 cm, assuming an 
initial separation <~ (Az),-//V). A lower bound on the Lorentz 
factor in this region can be obtained from the requirement that 
the optical depth to photon-photon pair production be less than 

1: 7 > 10 3 e ( 1/ 1 4 7V2 /4 (Afi)- 1 / 4 (l+z r ) 1 / 4 i/ 65 1/2 , where e M GeV is 

the energy of a test photon, 657/65 km s" 1 Mpc" 1 is the Hub- 
ble constant, and z r is the redshift (Woods & Loeb 1995). The 
solution exhibited in Figure 1 demonstrates that magnetically 
accelerated flows can attain the requisite high values of 7. 

Although our analysis can be generalized to arbitrary mag- 
netic field configurations, we have concentrated on the self- 
similar form for which we have exact solutions of the MHD 
equations (see VK). Real GRB outflows evidently have a fi- 
nite opening half-angle (estimated to be at least 5° on average; 
Meszaros et al. 1999), and thus are not accurately represented 
by our asymptotically cylindrical solutions. This discrepancy 
should not, however, affect our basic conclusions about the 
robustness of the magnetic acceleration mechanism for these 
flows. It remains a challenge for future work to find exact solu- 
tions with conical asymptotics. We note in this connection that, 
although the outflowing "shells" conceivably do not fill the en- 
tire solid angle into which they are ejected, the internal shock 
mechanism requires that multiple shells be ejected along any 
given direction (e.g., Kumar & Piran 2000). The guiding prop- 
erty of a magnetic field in an MHD-driven outflow provides a 
natural physical basis for this picture. 

In conclusion, we have shown that ordered magnetic fields 
can transform up to <~ 50% of the energy deposited by the cen- 
tral source of a GRB into kinetic energy of a collimated flow of 
baryons with 7 <~ 10 3 . This energy, in turn, may be converted 
into radiation by internal shocks. 

We thank J. Granot for helpful conversations. This work was 
supported in part by NASA grant NAG 5-9063 and by DOE 
under Grant No. B341495. 
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